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quarks: B ^ J/'ip {r]c) D and J (?7c) D tt 
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Abstract 

We evaluate the form factors describing the semileptonic decays J/fj [gf) D'^ I~nii 

within the framework of a QCD relativistic potential model. This decay is complemen¬ 
tary to JIf} ( 77 c) D^ 7r“ in a phase space region where a pion factors out. We 

estimate the branching ratio for these semileptonic and nonleptonic channels, finding 
BTi{^ Jjfj (gc) D+ £ vf) ~ Jjfj D+ tt”) = 3.1 x 10"® and 

Bn(^ gc vr-) = 3.5 x 10"®. 


PACS Numbers: 13.20.He, 12.39.Ki, 12.39.Pn,12.39.Mk 


In this article we study the meson decays 

( 1 ) 
( 2 ) 


50 ^ jj^ 

—>■ J/ifj (rjc) D~^Tr~ . 


These decays which may be classihed as “very rare”, at least in the Standard Model with 
mesons built solely of quark-antiquark pairs, should be searched for in present and future pro¬ 
grams at 5-factories. As a possible motivation consider preliminary studies of the inclusive 
B —> J/iIjX spectrum, indicating a slow J/ip hump [^, Q, which kinematically corresponds 
to an invariant mass mx — 2 GeV. Some hypotheses have been already suggested in order to 
account for such a phenomenon: in [Q the B^ ^ J/ijj A p decay is followed by the resonant 
A — p bound state, whereas in a possible explanation is the intrinsic charm content of the 
5-meson. In the latter case the decay proceeds through the B^{bd{cc)slow) —^ J/'^siowD*^ 
channel and it accounts for a BTZ ~ 10“^ provided that the intrinsic charm content fraction 
in the incoming 5-meson, is at least 1%. 

To corroborate the hypotheses in [§, 0, it is worth estimating the mechanisms for these 
decays in the framework of conventional heavy mesons picture as precisely as we can. In 
the 5° —*• J/pj {rjc) 5 ^° decays have been calculated in perturbative QCD; the branching 
ratios for these decays are estimated around and, therefore, too small to account 

for the Belle and CLEO data. Moreover in the possibility of production of a hybrid sdg 
meson with mass around 2 GeV is briefly discussed and, although the calculation of such a 
decay is difficult, a decay rate 10^ ~ lO'^ larger than the conventional mechanism for BTZ is 
expected. 

The nonresonant decay mode 5° —> J/pj {j]P) 7r~ would be interesting to analyze 

in this context, as it might provide a signihcant background to the decay process 5° — 
J/p) (pc) followed by 5*° — D'^ tt~ with a slightly off-shell 5*°. While a calculation 
from Erst principles is not available at the moment, a useful approximation might be the 
factorization approximation P] and, within this framework, the decay modes (2) would 
provide the crucial hadronic matrix elements needed to compute the relevant amplitudes. 


From a theoretical standpoint, semileptonic 5-meson decays with two hadrons in the 
hnal state represent a formidable challenge as they involve hadronic matrix elements of weak 
currents with three hadrons. They can be studied by pole diagrams, which amounts to a 
simplihcation because only two hadrons are involved in the hadronic matrix elements. This 
is the approach followed in some papers where these decays have been examined in the 
framework of the chiral perturbation theories for heavy meson decays Bi- This method 
has been successfully applied to decays with light mesons in the hnal state and it is based on 
an effective theory implementing both heavy-quark and chiral symmetry H, P, 1^ O . This 
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method allows to achieve, for systems comprising of both heavy (Q) and light (q) quarks, 
rigorous results in the combined mg —> oo, niq —0 limit. However the range of validity 
of this approach is limited by the requirement of soft pion momenta and it has been never 
applied, to our knowledge, to a hnal state with three charm quarks. The aim of this article 
is to examine the decays (|I|) in the framework of a QCD relativistic potential model [^] and 
to extend the kinematical range where theoretical predictions are possible. We shall present 
a detailed analysis of the form factors relevant for (2). Subsequently, the decays (@) will be 
considered. We do not include hnal state interactions in our calculation since no consistent 
way to compute them is presently available. It is clear that they can modify our numerical 


results ||r3[| . 

In three recent papers [|1^ an analysis of some semileptonic and rare H-meson 

decays into one and two light hadrons has been presented; it employs the relativistic potential 
model in an approximation that renders the calculations simpler. We wish to exploit here 
this approximation in the study of the B ^ J/'ip {rjc) D iu decay. 

Let us start with a description of the model (for more details see [| 
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In this 

approach the mesons are described as bound states of constituent quarks and antiquarks 
tied by an instantaneous potential H(r), which has a conhning linear behaviour at large 
interquark distances r and a Coulombic behaviour ~ —as{r)/r at small distances, with as{r) 
the running strong coupling constant (the Richardson’s potential is used to interpolate 
between the two regions). Due to the nature of the interquark forces, the light quarks are 


relativistic; for this reason one employs for the meson wave function T the Salpeter [1^ 
equation embodying the relativistic kinematics: 


+ R(r) 


vl>(f^ = M\l/(f) 


( 3 ) 


where the index 1 refers to the heavy quark and the index 2 to the light antiquark; M is 
the heavy meson mass that is obtained by htting the various parameters of the model, in 
particular the heavy quark masses which are htted to the values rric = 1452 MeV, mi, = 4890 
MeV, and the light quark masses ~ md = 38 MeV, = 115 MeV. The heavy meson 
wave function 4/jif(r) in its rest frame is obtained by solving eq. (|^) {H = heavy bound 
state); a useful representation in Fourier momentum space was obtained in and is as 
follows 

'0(/c) = 47rV Ma^ , (4) 

with a = 2.4 (1.6) GeV“^ for B, D (J/t/', rjc) mesons and k = \k\ the quark momentum 
in the heavy meson rest frame; this is the hrst approximation introduced in |^. At this 
point we would like to comment that the spectrum obtained depends only weakly on the 
light quarks masses. We will exploit this fact later when we employ = Aqcd- 
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The constituent quark picture used in the model is rather crude. There are no propagating 
gluons in the instantaneous approximation i.e., the Coulombic interaction is assumed to be 
static. Moreover, the complex structure of the hadronic vacuum is simplihed: conhnement 
can be introduced by the linearly rising potential at large distances, but the chiral symmetry 
and the Nambu-Goldstone boson nature of the tt’s cannot be implemented by the constituent 
quark picture. For these reasons, while there are good reasons to believe that eq. (|]) may 
describe the quark distribution inside the heavy meson, one cannot pretend to apply it to 
light mesons. Therefore pion couplings to the quark degrees of freedom are described by 
effective vertices (for more details see 0, 0,0). 

To evaluate the amplitude for semileptonic decays, it is useful to follow some simple rules, 
similar to the Feynman rules by which the amplitudes are computed in perturbative held 
theory. The setting of these rules is the main innovation introduced in |]^ as compared 
to [^. For the decays (p we draw a quark-meson diagram as in £g. and evaluate it 
according to the following rules: 

1) For the heavy meson H in the initial state one introduces the matrix: 




rriqmQ 


ii + "iQ p - #2 + rriq 


’mqUiQ + qi- q2 ‘2,mQ 2mq 


( 5 ) 


where mg and niq are the heavy and light quark masses, q^, q^ their 4—momenta and 
F = —* 75 , (/) for a. = 0“ (1“) heavy meson. The normalization factor corresponds to 


the normalization < H\H >= 2 mn and 


d^k 

W) 


\'4’H{k)\‘^ = 2mH already embodied in 


One assumes that the 4—momentum is conserved at the vertex HqQ, i.e. qi + q 2 = = 

hf-meson 4—momentum. Therefore gf = {Eg, k), q2 = {Eg, —k) and 


Eg + Eq = ruH . 


2) For the heavy meson H in the dual state one introduces the matrix: 

- 7"w7” . 


( 6 ) 

( 7 ) 


where H is as dehned in eq. (|^). 

3) To take into account the off-shell effects due to the quarks interacting in the meson, one 
introduces running quark mass mg{k), to enforce the condition 

E = \Jm^{k) + |fc|2 (8) 


for the constituent quarksQ. 

*By this choice, the average < rngik) > does not differ significantly from the value mg fitted from the 
spectrum, see [|^ for details. 
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4) The condition rriq > 0 implies the constraint 

0 < /c < /c, 


( 9 ) 


on the integration over the loop momentum k, where k^ax actually depends on the kinematics 
of the proces^ 

r 

( 10 ) 


J (2ir)3 

5) For the weak hadronic current one puts the factor 

7V,7V^.7^(1-75). 

The normalization factor Nq is as follows: 


( 11 ) 


m 


Nq=i 


(if q = constituent quark) 
Eg 


( 12 ) 


1 (otherwise) . 

6) Finally the amplitude must contain a colour factor of 3 and a trace over Dirac matrices. 

This set of rules can now be applied to the evaluation of the hadronic matrix element for 
the decays (|I]), corresponding to the diagram in £g. |T]; the result is: 



Figure 1: Feynman diagram for J semileptonic decay. J = J/ijj ,?7c- 


A>^ =< D+{pd) J\cY{l - 75)b\B%pB) >= 


(Pk 


(27r) 


U i}B{k) 'tp*j{k) i}*j^{k) 


0[kmax -k]Tr 


<b +md^b- + me'. 

2 rUb 2 nid 


4 +m. 


2 rriri 


2 rrir 


7 '‘(l -75) 


(13) 


tFor the processes induced by the b ^ u current, for instance, kmax — ms (see H, 1^ Q). Here 


kmax — 
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where T (—* 75 ) and J = J/'0(pj,e*) ( rjdpri) ) for the outgoing J/-^ {r]c) meson. Note 
the appearance of me' 7 ^ rric resulting from (|). The factor Af in ( p!3D accounts for the 
normalization of the hadronic current and of the heavy mesons: 


M = 


rricmb 


mdUib 


rndTUc' 


Eb{Eb - g°) Y EdTUB + md{mb - rud) V md{md - rric') + Ed{Ej -1713 + g°) 


mcrric' 


m% + — 2mcmc' + 2m^ + mjy — rrij — 2 {rnsq^ — Ed{Ej — 2mB + 2 g°)) 


(14) 


We introduce the various form factors for the B ^ J/'ijj D semileptonic decay: 


< D^iPn) J/'0(pj,e*)|c7^(l-75)6|50 (pb) >= * el 

{ + Ca,-) + q>‘{Ctp'h + C.q-') + p^,{C,p'i + Cjq’')] + 

+ i''"“^{D-iPj„qj3 + DlPJaPnq + Osq^pg/,) + -^(‘^“'’pBaPjaqqiDiPB + Diq’') + 

ruB 


_}_^y<^al3pB^pj^qp[DQp% + DjP^j + D^q^) \ , (15) 

mB J 

where q = Pb — Pd — Pj- Following we introduce also the form factors for the B ^ pc D 
semileptonic decay as follows: 


< D^iPo) Vc{Pv)\n'"{^ - 75 )&| 5 °(Ps) >= 

= i w+ {pd + pX + * {pr^-pnT + ir q>^ + 2 h PaP^qPDS ■ (16) 


It is useful to introduce the following variables: 


that satisfy 


s = {PD+Pjf 

t = ip- Pjf 

u = (p-Pd)^ , 

s + t + u = q‘^ + m^B + en^D + 


(17) 


The form factors tCi, r, h, Q, (i=l,...,7; j=l,..., 8 ) are functions of three independent 
variables. One can choose as independent variables s, g^, t or, alternatively, s, g^, Ej, where 
Ej is the J/-^ (? 7 c) energy in the i?-meson rest frame. The relations between the two sets of 
invariants are: 

t = m^B + ~ 2mB Ej 

g2 = s + m% - 2 mB {Ed + Ej) . (18) 
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The kinematical range is as follows: 


VI 

Oil Q 

t 

< {niB - mjf 

0 < 


< (jub -mj - y/i) 

^min — 

s 

<C. s 

_ ^max ? 


( 19 ) 


where 


^ min! max 

\{x,y,z) 



+ 11 ^ + — 2{xy + yz + xz) . 



( 20 ) 


From eq. (pTSf ) one can extract the different form factors by multiplying by appropriate 
momenta (see the Appendix for explicit expressions of all the form factors). 


The calculation of the trace in (^) is straightforward and is similar to those performed 
in 00 0 for similar processes. The evaluation of the integral is even simpler, because, 
although the three-body decay kinematics is rather involved, all the quarks in the loop 
of hg. g are constituent and there is no light quark propagating. The integration can 
be performed numerically but, unlike the semileptonic decays with two pious in the hnal 
state El. here the loop integration domain is not genuinely three-dimensional, due to the 
smallness of phase space. In fact, the calculation becomes simpler by inserting the wave 
functions {H = D, J/iIj, rjc) as they are in the Ff-meson rest frame in the relevant 

formulae, which is an approximation we perform and is justihed by the small value of the 
outgoing mesons’ 4-momenta in the processes (|l|). This approximation has been already 
incorporated into the formulae of the Appendix []. 


An important point to be stressed is the kinematical range in which the predictions of the 
present model are reliable. We cannot pretend to extend our analysis to very small meson 
momenta for the following reasons: hrst, as discussed in 0 0, when the outgoing meson 
momenta are small, the results of the model become strongly dependent on a numerical 
input of our calculation, i.e. the value of the light quark mass rrid- The numerical value of 
nid cannot be hxed adequately because the values of the quark masses were htted from the 
heavy meson spectrum, which is not very sensitive to the light quark mass (for more details 
see [0). Therefore the value of rrid has a theoretical uncertainty and the parameter rrid is 
the main source of error for the present calculation. We set the value of this parameter to 
the scale Aqcd0- 


^The form factors in (15 do not differ significantly from their actual value in the kinematical range 
( |l^ and within this approximation. 

Tn 10 the low-lying p-resonance provides the model with a cut-off at small s in the B —f tttt semileptonic 
decay form factors; here such a natural hadronic scale is absent. 
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Moreover the role of pole diagrams such as those studied in becomes relevant. These 
diagrams cannot be accounted for by the present scheme, which at most can be used to 
model a continuum of states, according to the quark-hadron duality ideas. The low-lying 
resonances, such as those studied in should be added separately^]. We expect a large 
contribution from the D* (see the discussion in [^). It is worthwhile to stress that in the 
present model the resonant production of a pion occurs through the long distance contri¬ 
bution depicted in hg. with T)h)+ mesons as intermediate states; all the diagrams are 
calculable in the Heavy Meson Chiral Lagrangean [j^j and they are found to be of the same 



J 


Figure 2: Feynman diagram for J nonleptonic decay in the QCD relativistic 

potential model. J = J/if^ ,r]c- 

In principle the partial widths F(i?o —> J iui) can be used to extract the relevant \Vcb\ 
CKM matrix element: due to the smallness of the phase space, the theoretical predictions 
are not expected to strongly depend on the specihc model employed to achieve the hnal 
results. We hnd BTZ{B^ —>■ J/'ip {rjc) i ug) ~ 10“^^, which is of course unmeasurable in 
the foreseeable future. 

Instead, we calculate the relevant formulae of the B^ ^ J ti~ nonleptonic decays. 
This channel is a background to —> J D*^ followed by the (almost on-shell) decay > 
7r“. The relevant amplitudes follow from eqs. ( 0 ) and (|T^). Numerically we get (for 
rud = 300 MeV ~ Aqcd): 

{ 1.98 X 10“® transverse polarization 

1.09 X 10“® longitudinal polarization , 

Bn(W r]c D+ TT") = 3.54 x 10“® . (21) 

^This is the reason why in the B* pole of the B ^ tt form factor is not reproduced in the |Fir| ^ 0 
region. The same remark holds for | ]l6[ about the resonances encountered at small s in the B —s- tttt 
semileptonic decay form factors, such as the p-resonance. 
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It is also interesting to compnte the differential branching ratios: 


dBU 

d\pj\ 


TBofi\Vciy:,?Gi \pj\ 
2567r3m^ Ej 


PSmax 

/ ds\qA{t,q^,s)\^ , 

^ ^min 


(22) 


where q^ = ml, f.^ = 132 MeV rgo = 1.6 ps, Vcb = 0.040, 14d = 1 — A = 0.22, ms = 
5.279 GeV, Gp is the Fermi constant and Sminimax are as in (^0[) . is the relevant amplitnde 
for the nonleptonic decay —>■ J/ip [pc) D~^ vr” given in eq. ([T3|) . The differential branching 
ratios for the nonleptonic decays we have stndied are plotted in figs. ^ and |^. 



Fignre 3: Differential branching ratio, in GeV for pc tt . The momentnm is in 

GeV. 


J/H^ Di££«reiDti«il Br^&obi&g Boibios 



Figure 4: Differential branching ratio, in GeV for B^ —> J/T D+ tt . The momentum is 
in GeV. 








We can therefore conclude that from an experimental point of view the i?—meson non- 
leptonic decay channel with three charm quarks in the hnal state represents a (barely) 
interesting process. We have investigated the —*■ J/'ip [rjc) D~^ i vi semileptonic de¬ 
cays and the —>■ J/if) {rjc) D~^ tt~ nonresonant nonleptonic decay channels by using 

the factorization approximation and the —*• J/ip {rjc) semileptonic decay form fac¬ 

tors. The calculation has been performed in the framework of a QCD relativistic potential 
model. The branching ratios for these decays are: BTZ{B^ J/xjj {rjc) D~^ £ ~ 10“^^, 
Bn{W ^ D+ IT-) = 3.07 X 10-® and Bn{W t]^ D+ tt") = 3.54 x 10“^ These 
theoretical hndings provide us with an order of magnitude for those processes and, in that 
respect, they do not seem to account for the slow J/xjj hump as indicated by the prelim¬ 
inary results of CLEO and Belle collaborations. Since the charmonium spectrum will be 
extensively studied at the 5-factories in the near future, it is important to conhrm whether 
the slow J/'^ hump exists. In that respect, a rehned measurement is needed. If the hump 
persists, it will be hard to hnd a consistent explanation within the conventional models. 
Thus a new scenario, like those discussed in i li, could be applicable. 
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in part by the Fund for the Promotion of Research at the Technion. 
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Appendix 

From eq. (PTSD one can extract the different form factors by multiplying by appropriate 
momenta. One getsll; 

/-fcmax dk k'^ , i^Bik) ij*j{k) Ipnik) r . ^ n 

Ci= U ^ ^^- md{md - me') + Ed{Ej - ms + q 

Jo /TT^ SVbmrimimrmtimh 


0£l I Ld' "‘C" ‘‘c'^b 
.3 r,™2 ^0 


\2EjmB{mB + q^) + Ed -m| - 2m|g° + {mb - md){mc - mc')q°+ 
Ej (-(ruft - md){mc + md) + mB{mB + g°)) + 

mB (mb{mc — me') + mcme' + 2md{mb — me + m^) — 3m^ — q'^'j 


+ 


rrid 


m\{me - me' + md) + EjmB{mb - me - 2md) + 


md{md + me){md - m^) + mBq^{me - me' + 2md) + 

Ej{mb - md)q^ + mdq^ - mb {{md - me'){md + me) + mB{mB + 2g°) + } , 


Oo = -2 


fkma. dk k'^ 'lpB{k) ij*j{k) 'lp*j^{k) p ^ 

/ A/ ^ ^ 5-5- md{md - m^) + Ed{Ej - mB + q 

Jo /TT^ bmi^m^memz,mb 


''B i"‘d"‘'c"‘‘e'^^b 
,2 


{Ed - mB) \mBmd{mb -me- 2md) + Ed{mj^ - {mb - md){me + md)) 

rkmax dk k‘^ , iJB{k) iJ*j{k) '0n(^) r / \ t-. / r-. n 

C3 = mB / „ 0 AT ^ —--r- md{md - m^) + Ed{Ej - mB + q 

Jo 6m%mymem^imb 

mBmd{2mb -me- m^ - Amd) + Ed{2m% - {mb - md){me + m^ + 2md)) , 


rkmax dk k'^ . 'lljB{k) '4’*j{k) ' 0 n(^) r / \ r-. /T-. On 

C4 = -mB / „ 2 a AT 2 —2-2- ^d{md - m^) + Ed{Ej - mB + q ) 

Jo /TT^ 8\'3m^^m^mcm^,mb 

[-2EdEbmB + mBmd{-2mb + me - me' + 2md) + Ed{{mb - md){me - m^ + 2md))] , 


rkmax dk k‘^ ^ i^B{k) 1p*j{k) ij*r,{k) v , . ^ n 

O 5 = 2m% / —— M — --- md{md - m^) + Ed{Ej - mB + q 

Jo Zn^ 8y 6m%m-^mcm‘^,mb 

[EdmB + md{mb - md)] , 

Cj = C, = -^C, , 

2 dk k^ ^ i>B{k) i>j{k) 'llj*D{k) \ , \ n / n 0^ 

Di = ml Af —---- md{md - m^) + Ed{Ej -mB + g°) 

JO /TT^ 8y Sm^^m'^mcme'mb 

[EdmB + md{mb - md)] , 


= b {j = 4,..., 8) is a consequence of the HqQ {H = B, D, J/b, Vc) couplings introduced in our 
model and of the nature of the quarks involved (all constituents): there is no way to generate higher powers 
in the meson 4-momenta in the loop of fig. 
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l-kma. dkP 'IpBi.k) 'lp*j{k) 'llj*j^{k) p 


Tn^rn^rricrnpmb 
mBmd{mb - rric - 2md) + Ed (m| - {mb - md){mc + md)) 


md{md - m^) + Ed{Ej - mB + g°) 


n _ _ dk k'^ 'il}B{k)'tp*j{k) i)*B{k) \ 

[{me - me') {Ed{mb - md) + mBmd)] , 


md{md - me') + Ed{Ej - mB + g°) 


-D 4 — Db — Dq — Dj — Dg — 0 , 


yj = i>B{k) iljj{k) iJD{k) 


md{md - me') + Ed{Ej - mB + g°) 


2 Jo 27 r 2 8V3ml m'^mcml,mb 

|- 4 Ej(mb - md){mB - q°) + mBmd 

[m| - 2me'{me + 2md) + md{3me + md) + mb{-me + 2me' + m^) + 2{Ejq° + g^) 
Ed -m\{me - 2me' + Amd) - 2md{-meme' + memd + ml + Ejq° + g^) + 
ImBmdq^ + mb (m% — Sm^g® + 2{—meme' + memd + ml + Ejq^ + g^)^ | , 


+ 


m\ rkmax dk k"^ . 'ipB{k) 'J>j{k) '0d(^) ' 

2 Jo 27r2 Sa/St^^ ^ 2 ^ ^ 2 , 


md{md - me') + Ed{Ej - mB + g°) 


mBm^mcm^'mb 

{mb - md) (-Ed(g° - mB) + md{md + me)) + m^ (md{q° - mB) + Ed{me + m^)) 


md{md - me') + Ed{Ej -mB + g°) 


rkmax dk k"^ ^ i^B{k) 'J>j{k) ' 0 d (^) ■ 

^ Jo 27r2 8y/3m\mlmcml,mb 

{EdfUb + Ebmd) mBq° + {mB - g°) (-Ej + 2Ed - mB) - q^ + {me + md) {me' - md) 


_ TT^ f>^rnax dk P ^ JjB{k) Jjj{k) ^^^{k) 

2 Jo 27r2 Sy/Sm'^ ^ 2 ^ ^ 2 , 


im^m^mempmb 


md{md - me') + Ed{Ej - mB + g°) 


[Ed{mb - md) + mBmd] 
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